Local dynamics and gravitational collapse of a self-gravitating magnetized Fermi gas. 
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We use the Bianchi-I spacetime to study the local dynamics of a magnetized self-gravitating 
Fermi gas. The set of Einstein-Maxwell field equations for this gas becomes a dynamical system in 
a 4-dimensional phase space. We consider a qualitative study and examine numeric solutions for 
the degenerate zero temperature case. All dynamic quantities exhibit similar qualitative behavior 
in the 3-dimensional sections of the phase space, with all trajectories reaching a stable attractor 
whenever the initial expansion scalar Ho is negative. If Ho is positive the trajectories end up in 
a curvature singularity that can be, depending on initial conditions, isotropic or anisotropic. In 
particular, if the initial magnetic field intensity is sufficiently large the collapsing singularity will 
always be anisotropic and pointing in the same direction of the field. 
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framework noted the possibility in which the pressure p^ 
parallel to the magnetic field overtakes the pressure p± 
perpendicular to the field, suggesting a magnetic collapse 
in this direction. This brings us to think that a general 
relativistic treatment of magnetized configurations will 
yield a singularity associated to the collapse having an 
extended, anisotropic, form along the direction of the 
magnetic field. Such type of singularities are denoted 
as "line" or "cigar" singularities, as opposed to isotropic 
"point -like" singularities. 

In order to investigate the local dynamics and col- 
lapse of a magnetized Fermi gas under General Relativity, 
we consider one of the simplest geometric configurations 
compatible with the anisotropic pressure that character- 
izes these sources: the Bianchi I model described in terms 
of a Kasner metric. The present paper generalizes previ- 
ous work along these lines [Hj]. 

The main justification of dealing with the simplified 
Bianchi I geometry is that the latter could provide a 
rough description of the local dynamics of a fluid ele- 
ment of the magnetized gas in the realistic configuration. 
If we consider local fluid elements far from the boundary 
of the configuration, so that this volume exchanges par- 
ticles and energy with the rest of the system (seen as a 
reservoir), then we can consider it like a local volume of 
a great canonical distribution associated with the whole 
gas. While it is evident that a lot of valuable information 
is lost by making such simplifications, we can still get a 
rough description of the effects of strong gravity on local 
physics. 

Once we consider the Bianchi I spacetime with the Kas- 
ner metric as the metric field associated with the mag- 
netized gas, the Einstein-Maxwell system of field equa- 
tions reduce to an autonomous system of four ordinary 
differential equations, leading to a 4-dimensional phase 
space that can be studied from a qualitative and nu- 

merical point of view using standard dynamical systems 

techniques. The physical and geometric dimensionless 
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I. INTRODUCTION 

White dwarfs and neutrons stars, as very dense objects 
found in great abundance, are astrophysical laboratories 
to test physical theories under intense gravity. In particu- 
lar, we can consider as a theoretical model the possibility 
of Magnetic White Dwarfs that could be endowed with 
extremely large magnetic fields, which, in principle, can 
be stronger than those measured in earthbound labora- 
tories. It is very interesting to study the behavior of the 
magnetic field lines, in its interplay with strong gravity, 
as such a magnetic star undergoes gravitational collapse. 
Various qualitative arguments show that under gravita- 
tional collapse the field lines squeeze together making the 
field stronger, even resisting collapse. 

Several models describing magnetic white dwarfs have 
been constructed (see [H, H, Hj] ) using Landau's well known 
argument setting appropriate mass limits for these ob- 
jects. Ideally, a consistent model of a magnetic white 
dwarf would require considering numeric hydrodynami- 
cal modeling of Einstein-Maxwell equations with (at the 
very least) axially symmetric configurations that would 
have to comply with the appropriate boundary conditions 
and equations of state of a degenerate magnetized Fermi 
gas. However, we can still obtain important information 
on the local dynamics of the self-gravitating magnetized 
Fermi gas by considering the evolution of such a source 
in a much more simplified spacetime geometry. 

The type of equation of state that we consider for 
a magnetized white dwarf is similar to that discussed 
in previous papers dealing with a magnetized electron 
gas [1, [H, @ (1968) and later developed for more gen- 
eral sources Q SHE, El (2000). These articles discuss 
equations of state for strongly magnetized gases endowed 
with anisotropic pressure. Previous work in a Newtonian 
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component of the shear tensor S, the expansion scalar 
7i, the normalized dimensionless magnetic field and the 
chemical potential, (3 and /x. All other quantities can be 
expressed in terms of these four basic quantities. The 
numerical examination of the system sheds light on the 
type of collapse singularities and their relation to specific 
initial conditions. 

The equation of state that we are using is strictly valid 
for densities of the order of (at least) 10 7 g cm -3 and up 
to 10 15 g cm -3 expected in compact objects from white 
dwarves to neutron stars. A gas of strongly magnetized 
and highly degenerate fermions in a neutron star is in a 
state that closely resembles superfluidity with near infi- 
nite conductivity (see page 291 [l3(). In these conditions 
the role of viscosity is minor, though one can still con- 
sider the possibility of dissipative or transport phenom- 
ena, such as dissipation of rotational energy in electro- 
magnetic and gravitational waves (see HEtll El). 

However, even if viscosity is not significant (at least for 
neutron stars), the main reason why we are neglecting it 
(and other dissipative effects) is to keep a mathemati- 
cally tractable problem. We feel that treating the case 
of thermal equilibrium is sufficient for a first approach, 
leaving the study of dissipative transport phenomena for 
a future work. 

The paper is organized as follows. In section II we 
present and discuss the set of Einstein-Maxwell equa- 
tions, the source of anisotropy and the most appropriate 
form of the equations of state for the magnetized Fermi 
gas. The qualitative dynamical analysis is carried on 
in Section-Ill, defining a set of dimensionless normalized 
variables, leading to a self-consistent autonomous system 
of four ordinary differential equations. The numerical 
and qualitative analysis of this system and the classifica- 
tion of the types of collapse singularities are dealt with 
in Section-IV, while the conclusion is given in section V. 
The main result from the numerical analysis is that once 
we allow for a relativistic strong gravity treatment it is 
always possible, for sufficiently large magnetic field ini- 
tial intensity, to obtain the anisotropic "cigar" type of 
collapse singularity as hypothesized in previous work Q 
carried on along an intuitive Newtonian framework work. 



II. KASNER METRIC WITH ANISOTROPIC 
PRESSURE 



The Kasner metric is among the simplest metrics com- 
patible with the anisotropic pressure associated with a 
magnetized source. This metric is given by: 



O and the shear tensor a% take the forms: 



ds 2 



-c 2 dt 2 + A 2 (t) dx 2 + B 2 {t) dy 2 + C 2 (t) dz 2 . (1) 



It is associated with a "non-tilted" Bianchi-I space time 
[III. For a comoving 4- velocity u a = Sf, where u a u a = 
— 1, the 4-acceleration vanishes and the expansion scalar 
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We consider as the source for this metric the following 
stress-energy tensor: 



T b a = (U + P)u a u b 
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n 



6> 



P = P 5—- (4) 



where 11^ is the anisotropic pressures tensor, U the en- 
ergy density, P the pressure, B the magnetic field and 
A4 the magnetization, all them are functions of the time. 
Notice that the anisotropy is produced by the magnetic 
field B. If this field vanishes, ie: B = 0, the stress-tensor 
reduces to that of a perfect fluid tensor with isotropic 
pressure P = p. In the general case B ^ the tensor 11^ 
has the form 



n£ = diag[n,n,-2n,o], n 



BM 



K = 0- (5) 



The Einstein field equations (EFE) associated with the 
Kasner metric and the stress-energy tensor are 
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where a dot denotes the derivative with respect to the 
proper time of fundamental observers and k = 87rG/c 4 . 



From the balance equation T' 



a b 



and the Maxwell 



equations F ab . b = and F[ ab;c ] = 0, wc further have 



U + (p + U)Q - BM{ 
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Since we need to construct a self-consistent system of 
first order differential equations that can be solved nu- 
merically, it is convenient to eliminate first and sec- 
ond derivatives of the metric functions in the Einstein- 
Maxwell equations in terms of the expansion scalar and 
the components of the shear tensor. Proceeding along 
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these lines, we combine equations ([2]),(|3]), ([6]) and (J7]) 
to eliminate the functions A, B, C and their derivatives 
A, A, B, B, C, C. After some algebraic manipulation we 
arrive to the following constraint 



B 2 

- (E a ) 2 = (E 2 ) 2 - E a E z + — =kU, 

3 

plus the following set of 5 differential equations: 



U + (U+p)O-BM(-O-Z z )=0, 
tv = ~^BM-^ V Q, 

o 

e z = ^BM-z z e, 

3 e 2 f 
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(3 = -^(3E 2 -29), with (3 = B/B C 



(9d) 
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where E 2 = a z is the independent component of the 
shear tensor. While the shear tensor can be fully deter- 
mined by this single quantity, it will be convenient for 
our numerical calculations further ahead to use two com- 
ponents of this tensor. 

The functions U, p and M. are now given by the equa- 



tion of state for the gas; 

p = \T p (p,n), BM = \(3T M (p,n), U = \T v {f3,n). 

. ( 10 ) 

The r functions depend of the parameter (3, which is 
the magnetic field normalized by B c , and by the chem- 
ical potential, /z, normalized with the rest energy, both 
dimensionless quantities. The constants B c and A are: 



B c = 



m 2 c 3 mc 
A = 



4tt 2 A3' 



(11) 



while A c = h/mc is the Compton wavelength. If we 
consider an electrons gas, then m = m e and A = 
3.86.10- n cm. Nevertheless, B c = 4.414 x 10 13 G (for 
electrons too) which is a well known critical value of a 
magnetic field. In all neutron or white dwarf stars older 
than few second after formation, one can neglect the 
thermal contributions to the pressure and energy den- 
sity; thus we can set p(/3,/z, T) = p(/3,/z), the same for 
U and M.. Typical white dwarf temperatures satisfy 
kT << -©Fermi, where ©Fermi is the Fermi kinetic en- 
ergy. So, the thermal disorder kT is not responsible for 
the pressure, the energy density nor the magnetization 
in pop . It is now convenient to introduce the form of the 
r(/3. m) functions for the degenerate case (T = 0), these 
are Q: 
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while a n = fiy/ fi 2 — 1 — 2n/3, b n = In 
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where a n — 2 — 6o n , n = 0, l,(.12c) 
> 2 -l" 



s = I 



2(3 



(12d) 



where I[X] denotes the integer part of its argument X. 



III. DIMENSIONLESS VARIABLES 



Consider now the following variables: 

TT 8 d Id 
H 



3 dr Ho c dt 
and the dimensionless functions: 



(13) 



where S y and S z are related to the yy and zz components 
of the shear tensor, while f2 is related to the magnetic 
field. The new time r is a dimensionless time (or "loga- 
rithmic" time). The quantity H(t) (because of (|3J)) will 
have dimensions of cm" 1 and the sign of r becomes de- 
termined from the sign of Hit), However, we have chosen 
ft = /3 and k\ = 3ff 2 in order to avoid the presence of 
annoying constants in the system of equations. Inserting 
the equations of state (TIT))) and the new definitions (|13[) . 
([T4"| into fS]) and © we obtain the constraint: 
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plus the system 



-^r^ - 3S y n, 
2f3v M - is z n, 
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(16a) 
(16b) 

, (16c) 
(16d) 



[{2H - S Z )(T M + 2T Ui0 )P - 3H(T P + Trj)} 

(16e) 



Notice that, as opposed to cosmological sources and 
models [ijj where H n = 0.59 x 10 _28 cm _1 would play 
the role of the Hubble scale constant, in our magnetized 
Fermi gas we have Hq — 0.86 x 10 _12 cm _1 , which 
is a much smaller length scale. This is logical and 
consistent because it indicates that our simplified model 
is examined on local scales smaller than cosmic scales. 
The scale 1/H ~ 1.15 x 10 12 cm is the order of the 
distance of an astronomical unit. 

The results mentioned in references @, H, Q for the 
electron gas show that for an intense magnetic field, of 
the order of the critical field B c , all the electrons are in 
the ground state of the Landau level n = 0, and con- 
sequently we have p±_ =0. It is an interesting issue to 
study how the electron gas evolves in this case, in which 
the functions T(/3, p) are simplified considerably, taking 
the following forms: 
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a = fiyj )i 2 - 1, b = ln(fi + ajfi), c = 0, a = 1. 

Thus, substituting (fl~7|) into (fT6|) yields a self-consistent 
system of five ordinary differential equations (|16[) . with 
the unknown functions /3,H.,S y , S z and fi, and the 
constraint (|15p . which only admits a numerical solution. 
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FIG. 1: Numerical solutions for the magnetic field (3(t) and 
the expansion scalar H.{t). Here we got three different initial 
conditions S^O) = 0,S V {0) = -1,S Z (0) = l,/*(0) = 2, and 
/3i(0) = 10~ 5 ,/3 2 (0) = 5.10- 5 ,/3 3 (0) = 10~ 4 respectively, for 



the initial magnetic field. Then because of (|18[) . the initial 
expansions should take the values 7ii(0) = Ttzfi) = T~ls(0) = 
—4.82. The numerical results for the collapse times are ri = 
0.1059, r 2 = 0.1058, and r 3 = 0.1057 respectively. 



For non-tilted Bianchi-I models there is only one inde- 
pendent component of the shear tensor, which we are 
taking to be S z , however it is useful to use also the S v 
component for several numerical calculations. Also, the 
form of the components of the shear tensor determines 
the form of the metric coefficients. From <(2])— ([3]) we get: 



c 



= (S z +H). 

(20) 



where S x 



-S y - S z 



From the equations of state (fTU|) it follows that the 
chemical potential must satisfy fi > 1, which is typi- 
cally correct for systems with densities of the order ~ 
10 7 gm/cm 3 or larger. For white dwarfs or neutron stars 
the chemical potencial takes values around V3 ~ 1.732. 
Since U > 0, then the chemical potential fi > 1 and from 
(|15p we obtain the constraint equation: 



(s v ) 2 ~{s z ) 2 - s v s z + m 2 



3Ttj > 0. 



(18) 



so that our physical 5-dimcnsional phase space is re 
stricted by the relations: 

,2 > (s y ) 2 + (S z ) 
fj 2 > 1 + 2/3, 



m 1 



s y s z 
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P > 



3r r 
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(19b) 
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IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

Using and (TT3"|) we can express the local average 
volumcn V = ABC in terms of H and the dimensionlcss 
time r as 



V(r) = V(0) exp 3 



Hdr 



(21) 



clearly illustrating why the time r is known as a "log- 
arithmic" time. Notice that the sign of H(t) denotes 
expanding (> 0) or collapsing (< 0) local volumes. 
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A. Singularities 

Since the equations of state that we are considering are 
associated with compact objects of very high density (at 
least ~ 10 7 gm/cm ), the evolution range of the models 
for diluting lower densities is not physically interesting 
and will not be pursued any further. This means that 
we will only consider the collapsing phase of the models, 
so that we will only need to examine initial conditions in 
which the initial expansion Ho is negative. We tested nu- 
merically the models using a wide range of different initial 
conditions covering the full range of physically interesting 
values for compact objects, from white dwarf to neutron 
stars, for example: fiQ = 2, corresponding to densities 
of ~ 10 7 g/cm 3 , while (3q — 10~ 5 represents magnetics 
fields of 10 s G. Together with Ho < 0, we considered in 
particular: Sq = 0, ±1 and Sq = 0, ±1, corresponding to 
cases of zero initial deformation and initial deformation 
on the y or z directions. 

As long as Ho < the models exhibit a general col- 
lapsing behavior H — ► — oo, independent of the initial 
values of other functions (see an example in figure^)). 
In all collapsing configurations the magnetic field inten- 
sity diverges to infinity regardless of its initial value. In 
figure ([T]) we show several numeric curve solutions for dif- 
ferent values of the initial magnetic field. We can see 
how for increasing intensity (labeled by (31 < (32 < (33 
in the figure) the expansion shows a faster decay to — oo. 
Consequently, whenever we increase the initial magnetic 
field the collapse times decrease. 

We consider initial condition with initial shear defor- 
mation not necessarily in the direction of the magnetic 
field, so that anisotropic singularities of type "cigar" 
emerge along any of the three axes, that is, a collapse that 
can be parallel or perpendicular to the magnetic field. An 
examination of all these cases reveals that the collapse 
state strongly depends on the magnitude of the initial 
shear. For example, if we have a large initial shear (de- 
formation) in the x direction, say: Sq 3> Sq , Sq , (3q , /j,q , 
then a type cigar singularity emerges along the x direc- 
tion. This is shown in figure©, illustrating (by means 
of the set of equations (|2H)) ) that the metric function A 
tends to infinity, while the other metric functions, B and 
C rapidly fall to zero. In general, the initial configuration 
of the system and the initial values of the shear tensor 
determinate the privileged direction of the anisotropic 
collapse. 

The numerical trials also show that there is always 
a threshold value for the initial magnetic field intensity 
that influences the direction of the "cigar" type singu- 
larity. This is illustrated in figure©: if we increase the 
initial magnetic field intensity to (3q — 1 then we will ob- 
tain a "cigar" type singularity along the z direction for 
any (3q > 1. However, (3q = 1 ~ 10 13 G, which is not 
a physically realistic value in magnetic white dwarfs but 
could be reasonable in a primordial magnetized universe 
model. 

An isotropic "point" singularity can always emerge for 
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FIG. 2: Trajectories of the functions S l + Ti for i = x,y,z. 
We take the following initial conditions: S x (0) = 1, S y (0) = 
0, S*(0) = -1, n(0) = 2, /3(0) = 0, and H(0) = -3.43. This 
system shows an initial deformation in the x direction with 
the magnetic field pointing in the z direction. The collapse is 
in the form of a "line" singularity in the x direction. 



configurations with zero initial deformation, i.e. Sq = 
&o = = an d A) = 0. However, even with zero 
initial deformation there is always a threshold value of 
initial magnetic field intensity for which the singularity 
becomes anisotropic along the z direction. 



B. Phase Space and Critical Subspaces 



As mentioned previously, the Einstein-Maxwell sys- 
tem can be written as an autonomous system associ- 
ated with a 4-dimensional phase space in the variables 
(S z , (3, fi,H). Notice that T, y can always be found if we 
determine S z = X, which is the only independent shear 
component. Considering the four primordial functions 
the system will take the form: 



(U + p - -BM)S - BME, 



u 

E = -j;BM - - Br 

e 



K (BM + -(U-p))-Q 2 , 
^/3(3E-26). 



(22a) 
(22b) 
(22c) 
(22d) 
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where £ = a z . If we work with the dimensionless func- 
tions defined in ([14")) this system becomes: 

S* = 2(3T M - mS z , (23a) 

n, T = f3T M + ^{T v -T p )-?,H 2 , (23b) 
£ T = 2f3{S z -2H), (23c) 

M,r = -^—[(2n-s z )(r M + 2r u ^)p-m(r p + r u )]. 

(23d) 

Note that only change the equation for H(t), but we 
can reduce the equations (123bj) or (I16c|) to the constraint 
(fT5)l . therefore both systems are equivalent. Then, p6a[) 
is only necessary for the computation of the metric coef- 
ficients. 

In figure ((3|) we represent a 3-dimensional section of the 
phase space (S z , /3, /i) with different curves for several ini- 
tial conditions. As shown in the figure, the initial value 
of the expansion Hq determines the global evolution of 
the numerical solutions curves. All curves starting at 
t = with initial expansion Ho = — y kA/3 converge 
into the stable attractor marked as "a" , whereas if we 
choose Hq — \J kX/3 and start at r = 0, the curves evolve 
towards the anisotropic singularity. Setting the left hand 
sides of the equations in ([23| to zero and solving the alge- 
braic system we find the set of critical points associated 
with this system, including the stable attractor marked 
as "a" . These points are given by: 

a = {S z = 0,/3 = 0,fi = 1,H = 0}, (24) 

For the 4-dim phase space (S z , (3, n,H) we have 4 pos- 
sible 3-dim sections of the same phase space. We have 
computed this numerical solutions with similar results. 

V. CONCLUSION. 

We have presented a model based on the dynami- 
cal description of a local volume of a magnetized, self- 
gravitating, Fermi gas in the basic Landau level n = 0. 
Since we are considering complicated equations of the 
state, we have worked with the simplified form of the 
Einstein-Maxwell equations that follow by assuming a 
Bianchi-I space time represented by a Kasner metric, 
whose source of anisotropy is just the magnetic field. 
This simplified spacetimc geometry provides a convenient 
toy model for a rough understanding the local collapsing 
behavior of the type of matter found inside a magnetized 
star like a white dwarf or a neutron star. 

The relevance of the present paper emerges from 
our study of the collapsing singularities, which can be 
isotropic point-like or anisotropic of type cigar. Point 
singularities emerge under very special initial conditions 
of zero magnetic field, zero shear deformation or both. 
Cigar type singularities can also be obtained in all direc- 
tions, depending on the initial values of the shear defor- 
mation. However, for an initial magnetic field intensity 




FIG. 3: Trajectories in a section of the phase space (S z ,/3, fi) 
for four different initial conditions. The shaded fi — 1 plane 
is bounded by -4 < S z < 4 and < j3 < 0.2. The stable 
attractor is the point marked by a, and the curves 6, c, d and 
e are numerical solutions on the 3-dim section of the phase 
space. All the trajectories here fall from r = toward the 
stable attractor a. Similarly, from r = (empty spheres) all 
them scape toward the singularity. 



having a sufficiently large value the end singularity al- 
ways becomes of type cigar in the direction of the field. 
This result is important because the value of the mag- 
netic field determines the type of collapse and this is in 
agreement with the non-relativistic previous paper which 
examined the collapse of this type of magnetized gases 
within a Newtonian framework 

As discussed in [lj| by Collins & Ellis, orthogonal 
Bianchi models like the one we are considering arc glob- 
ally hyperbolic and only present a single singularity, so 
that hypersurfaces of constant time (orthogonal to the 
4- velocity) are global Cauchy hypersurfaces and every 
point in spacetime can be causally connected to the lat- 
ter. The 4- velocity is a geodesic field and the singularity 
is marked by a specific constant time value, so that it is 
non-timelike and every event in spacetime can be causally 
connected to this singularity (in particular by the time- 
like geodesies that are integral curves of the 4-velocity 
field). 

Now, in this article we are only considering a collaps- 
ing regime from an initial hypersurface of constant time. 
Thus, every future directed timelike curve (geodesic or 
not) starting at any initial Cauchy hypersurface of con- 
stant time will terminate in the collapse singularity. Un- 
der these conditions, this singularity is obviously cen- 
sored. 



7 



Acknowledgments The authors acknowledge the Office of External Activi- 

ties of ICTP for its support through NET-35. The au- 
We gratefully acknowledge very helpful through the fi- thors are also indebted to Dr. Hugo Perez Rojas by the 
nancial support from grant PAPIIT-DGAPA-IN117803. permanent interest in this work. 



[1] Christian Y. Cardall, Madappa Prakash and James 
M.Lattimer 2001 The Astrop. J ,554: 322-339. 

[2] A.Broderick, M.Prakash and J. M.Lattimer 2000 The As- 
trop. J, 537: 351-367. 

[3] I.Bednarek, A. Brzezina, R. Manka, M. Zastawny-Kubica 
2005 Theoretical model of a magnetic white dwarf 
Preprint astro-ph/ 0212483 vl. 

[4] V. Canuto and Hong-Yee Chiu 1968 Phys. Rev. Vol. 173, 
No-5, 1210. 

[5] V. Canuto and Hong-Yee Chiu 1968 Phys. Rev. Vol. 173, 
No-5, 1220. 

[6] V. Canuto and Hong-Yee Chiu 1968 Phys. Rev. Vol. 173, 
No-5, 1229. 

[7] Felip.R. Gonzalez, Cuesta.H. J. Mosquera,Martmez. A. Perez, 

Rojas.H. Perez 2005 Chin. J. Astron. Astrophys, 5 399, 

( Preprint |astro-ph/0207150[ ) . 
[8] Martinez. A. Perez, Rojas.H. Perez, Cuesta,H. J. Mosquera. 

2003 Europ. Phys. J. C 29 111-123. 
[9] M.Chaichian,S.S.Masood,C.Montonen,A. Perez Martinez 

and H.Perez Rojas 2000 Phys. Rev. Lett 84, 5261. 
[10] Somenath Chakrabarty 1991 Phys. Rev. D Vol. 43 No. 



2. 

[11] Somenath Chakrabarty 1996 Phys. Rev. D Vol. 54 No. 
2. 

[12] A.Ulacia Rey, A.Perez Martinez, Roberto. A. Sussman 
2007 Int. J. Mod. Phys. D, Vol. 16, Nos. 2-3, p 481-487 
(Preprint gr-qc/0605054 ). 

[13] Stuart L.Shapiro, Saul A. Teukolsky. 1983 "Black Holes, 
White Dwarfs, and Neutron Stars" ed John Wiley & 
Sons, Inc. 

[14] Stuart L. Shapiro. 2005 Preprint [gr-qc/050 9094 vl. 

[15] Mattew D. Duez, Yuk Tung Liu, Struart. L. Shapiro, 
and Brason C. Stephens. 2004 Phys. Rev. D Vol. 69, 
p.104030. 

[16] J. Wainwright, G. F.R.Ellis 1997 Dynamical Systems in 
Cosmology, ed Cambridge University Press. 

[17] Charles W. Misner, KipS. Thorne, John. Archibald 
Wheeler 1973 Gravitation, ed WH. Freeman and Com- 
pany, NY. 

[18] C.B.Collins, G.F.R. Ellis. 1979 Phys. Rep, Vol. 56, No 2, 
65-105. 



